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Abstract 

In this paper we study the effect of non-magnetic impurities in spin-Peierls 
system CuGeOs within the framework of non-linear sigma model plus topo- 
logical term. We show that local moments are induced in both the Zn-doped 
and Si-doped CuGeOs compounds. Effective low energy theories for the 
impurity-induced local moments are derived in both cases where interest- 
ing differences between the two cases and between the Zn-doped Heisenberg 
two-ladder system are pointed out. The low-energy physics of the Zn-doped 

Heisenberg three- leg ladder system is also discussed. 
PACS Numbers: 75.10.Jm, 75.10.-b, 75.40.Mg 
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Recently, there are lot of interests in the study of quasi- one dimensional spin 1/2 Heisen- 
berg spin systems, including the spin-Peierls system CuGeO^ |]J and Heisenberg ladder 
systems Sr n Cu n+ \02n+i • In particular, a lot of efforts have been spent on the study 
of impurity effects on CuGeOz and SrCu20s compounds where surprising results were ob- 
served. In CuGeOs, which is an inorganic spin-Peierls system, it was observed that local 
magnetic moments were induced when a few precent of Cu is replaced by Zn or when 
Ge is replaced by Si ||. In the case of Cu\- x Zn x GeOz, it was discovered that the spin- 
Peierls state collapses at around x = 0.03, and is replaced by a new magnetic state which 
was interpreted as a spin-glass-like state 0] or long-range antiferromagnetic state ||. In 
CuGei-ySiyO^, long-range antiferromagnetic order co-existing with spin-Peierls state was 
also observed with very small amount of disorder, when y = 0.007 f7]||. Similar 'disor- 
der induced' local magnetic moments were also observed in the two-leg ladder compound 
Sr{Cui- x Zn x )20^ when Cu ions are replaced by Zn ions || and long-range antiferromag- 
netic order was observed at low temperature at x ~ 0.07 |IIJ]. The possibility of co-existing 
spin-Peierls and long-range antiferromagnetic state and formation of local moments in Zn- 
doped two- leg ladder system has been investigated by bosonization technique |TI|jr2]] . It 



was also suggested that the low-energy physics of the Zn-doped Heisenberg ladder system 
can be described by a S = 1/2 random J Heisenberg spin chain where the nearest-neighbor 



coupling J has both random sign and magnitude |L3 . 

In this paper we shall discuss the effect of non-magnetic impurities in the spin-Peierls 
system within the framework of non-linear sigma model (NLaM) plus topological terms 

L = \d,nf + *A(|n| 2 - 1) - i^n^n x d u n), (1) 
g 8tt 

where n is a unit vector and g ~ 2/5 = 4. The model is believed to be the correct starting 
point for Heisenberg spin chains [|14],[nj. We shall show that the formation of magnetic 
moments in both Zn-doped and S'i-doped CuGeO^ compounds can be explained within 
the NLaM treatment. The corresponding low energy theories for finite concentration of 
impurities will be derived where similarities and differences between the two cases and 



between the Zn-doped Heisenberg two-ladder system will be pointed out. The case of 
Zn-doped Heisenberg three-leg ladder system will also be discussed. We shall argue that 
the experimentally observed 'disorder induced' magnetic states in different compounds are 
consistent with our effective theories. 

In the NLaM formulation, a dimerzied S = 1/2 spin chain with alternating interaction 
Ji t i+i = J[l + 7(— 1)*] can be described for small 7 by a NLaM with with topological term 
9 = 7r(l + 7) |15|| . For an infinite chain, there is a spin gap A g ~ J7 2 / 3 in the excitation 



spectrum |T(J. We shall assume that the only effect of substituting Cu by Zn is to remove 



spins randomly from the spin chain, leading to a collection of open spin chains with average 
length ~ x -1 , where x is the Zn concentration. 



The properties of open spin chains have been studied by the author [IT) where formation 
of boundary excitations (end states) and their topological characters were discussed. For 
open S = 1/2 dimerized spin chain with even number of sites % = 1, 2m, it was found that 
localized 5=1/2 excitations exist at the ends of spin chain when 7 > 0, but are absent 
otherwise. The existence of end states when 7 > can be understood most easily in the limit 
I7I — > 1. In this limit, the spin chain is completely dimerized and no end states appear if J^+i 
is finite for odd i and zero for even % (7 < 0). In the other case when J^+i is nonzero when i 
is even (7 > 0), the first and last sites are decoupled from the rest of the spin chain forming 
the 5 = 1/2 end states. For I7I small, the end spins are not completely decoupled from 
the spin chain and end states which decay into the spin chain with finite correlation length 
£ ~ A^ 1 ~ |7|~ 3 / 2 are formed fLTj . For the randomly doped Cu\^ x Zn x GeO^ compound, 
it is expected that open spin chains with both even and odd number of sites are formed. 
However, the polarization of the dimerization 7 is not determined a priori but is determined 
by the underlying dynamics of the broken spin chain. In particular, we expect that even 
chains will always dimerized with 7 < to maximize the number of short-range spin-singlet 
bonds and no end states are formed. However, the situation is different for spin chains with 
odd number of sites where the ground state is expected to be a spin doublet with total spin 
Stot = 1/2. To respect parity, we expect that the spin chain will be dimerized in a symmetric 



way leaving an 'extra' spin S = 1/2 as a localized state at the middle of spin chain. Notice 
that in this case, the spin chain has opposite polarization of dimerization at two ends of the 
spin chain, and the localized S = 1/2 state at the middle of the spin chain is just a soliton 
solution at the domain wall seperating the two regions of polarization ||18|| . 

The existence of localized S = 1/2 state between two regions with different polarization 
of dimerization can also be understood from the topological consideration of the NLaM. In 
this description, the spin chain is equivalent to a a-model with different topological angles 
6l = tt(1 — 7) and Or = tt(1 + 7) on two sides of the spin chain, where 7 > 0. The 
qualitative properties of this model can be undertsood using a large-N expansion in the 
CP/v-i representation of the NLaM, where spin excitations are represented by charged- 
bosons in a one- dimensional universe of scalar QED jl9) . 

L ~ — 
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\(d, + iA,)Zf + AJ|Z| 2 ] + ±F% ~ (2) 



where Z 



is a two component spinor boson field, = dnA u — d u A, t . e is the electric 



charge carried by the bosons and the sign of the charge is in fact a sub-lattice index |T7 
The #-term is given by 9{x) = &l(r) on left(right) sides of spin chain and is interpolating 
smoothly between the two regions. In this description, the #-term gives rise to a background 
electric field Eb(x) = (8(x)/27r)e, and the ground state is obtained by minimizing the energy 
of the system in the presence of £7&(#) [TIJ. The electric field energy of the system can be 
reduced by nucleating a boson of charge — e out from vacuum and localized at the boundary 
region where 7 changes sign. The total electric field is the sum of the background electric 
field and electric field radiated from the charge boson. The resulting electric field remains 
unchanged at the left-hand region, but is reduced to E = ((l + j)/2 — l)e = —(l — j)e/2 on 
the right-hand region. The resulting change in energy is of order A g — e r y(L/2), where A g 
is the energy cost in nucleating the boson, and L is the size of system. The energy of the 
system is always reduced for large enough L. The nucleated boson in this theory is precisely 
the localized 5 = 1/2 state at the middle of spin chain. 



The localized spin states between regions of spin chain with different topological angle 
9 can also be used to understand the local moments generated in the CuGei_ y Si y 3 com- 
pound. In this case, the spin chain is not broken into open chains by Si. However, replacing 
Ge by Si changes the coupling between spins close to the Si ion, or the elastic constant 
governing the lattice distortion, and may lead to a local broken symmetry where one polar- 
ization of dimerization is favoured over the other for the region close to the Si impurity. For 
finite concentration of Si randomly distributed in space, the favoured polarization is ran- 
domly distributed, and the system can be described as a dimerized spin chain with domain 
walls formed between 'mismatch' Si impurities, or equivalently, a NLaM with topological 
angle 9(x) = 7r(l ±7) randomly distributed in space. Localized S = 1/2 state will be formed 
at boundaries between regions with different 9 as demonstrated above. 

The appearance of localized magnetic moments leads to Curie behaviour in uniform 
susceptibility x ~ ( X ,y)/2T a t low temperature T < A g and low doping, for both 
Cui_ x Zn x Ge0 3 and CuGe^ySiyOs compounds. Notice that the concentration of mag- 
netic moment is on average, half of the concentration of impurities in both cases. At lower 
temperature and higher concentration of impurities, interaction between magnetic moments 
becomes dominant and new magnetic states may develope in the systems. In the follow- 
ing we shall derive the effective Hamiltonian for the impurity-induced magnetic moments. 
To derive the effective Hamiltonian, we divide our system into two parts, the bulk, dimer- 
ized system and the impurity-induced local moments, the bulk system is described by the 
Hamiltonian with the local magnetic moments removed, and the total Hamiltonian is 

H — H a + J' Si'.Sj, 

<i',j> 

where H Q is the Hamiltonian for the bulk system, SV are the localized moments, j's are the 
sites in the bulk system next to i' , and J' < J is some effective coupling between the local 
moments and spins in the bulk system. The effective interaction between the local moments 
can be obtained by integrating out the bulk system. To second order in J', we obtain an 
effective interaction between local moments coming from exchange of bulk spinwaves, 
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Heff ~ ^2 -Ji'j'Si'-Sj', (3a) 

i'j' 

—> — * 

where SV and Sj* represent the localized magnetic moments, and [T3| 



j.,. (J') 2 Y,<S t .S j >, (3b) 

ij 

where i's and j's are the sites next to the local moments i' and j', and < Si.Sj > is the 
zero-frequency spin-spin correlation function between sites i and j in the bulk system. In 
particular, for large distance = |r£ — fj-| and in the spin-Peierls state, we expect 



<^4->~(-l)^'e-^, (3c) 

where £ ~ A~ x is the correlation length of the bulk spin system, and < S^Sj > is positive 
when i and j are on the same sublattice, and is negative otherwise, because of the underlying 
antiferromagnetic correlation in the system. As a result, Jyji is negative between local 
moments on the same sublattice, and is positive otherwise. Notice that intra-chain coupling 
can also be included in the effective Hamiltonian by including intra-chain coupling in the 
bulk Hamiltonian H Q . The qualitative result (0) is not modified except that the correlation 
length £ is anisotropic with a much larger value along chain and smaller value between 
chains. 

Now we may apply Eqs. (ff) to the Cu\- x Zn x GeOz and CuGe\- y Si y O^ systems. For 
the Zn-doped system, inter-chain coupling between localized moments is extremely small 
because the spin chain is broken, and the induced local moments are located at different seg- 
ments of the broken chain. The ground state behaviour of the system is mainly determined 
by intra-chain coupling, resulting in an effective random (but unfrustrated) Heisenberg an- 
tiferromagnet model in 3D. The coupling is random because the locations of the local 
moments are random. However the system is unfrustrated because although the magnitude 
of the coupling J^'s are random, the sign of the effective interaction is not: local moments 
on the same sublattice always interact ferromagnetically whereas local moments on different 
sublattice always interact antiferromagnetically. The classical ground state of the impurity 
system is always the Neel state which is unfrustrated. 
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For the .Si-doped system, the consideration is similar except that inter-chain interaction 
between local moments exists and dominate the physics at an intermediate temperature 
regime when intra-chain coupling is not as important. Because of the topological nature of 
the domain walls, we find that the local moments are always arranged in a staggered way 
where the local moments adjacent to each other are always locate on opposite sub-lattice 
and interact antiferromagnetically, i.e. the impurity-induced local moments can be described 
by a model of ID antiferromagnetic spin chain with nearest neighbor interaction Jy which 
is always positive but with random magnitude. At lower temperature when intra-chain 
coupling becomes important, the system is described by an effective 3D random Heisenberg 
antiferromagnet as in Zn-doped case. 

The theoretically suggested similarity and difference between the Zn-doped and Si- 
doped CuGe0 3 compounds seem to be in qualitative agreement with what is observed 
experimentally. At very low temperature when intra-chain coupling becomes important, 
both systems are described by a model of 3D random Heisenberg antiferromagnet which, in 
the continuum limit, can be described by a (3+l)d NLaM with random spin-wave velocity 
c(x) and random spin stiffness p s (x), 



and has an ordered phase when the effective randomness is weak enough. This long-range 
ordered antiferromagnetic phase seem to be observed in both compounds A crucial 

difference between the two compounds is that in the Zn-doped compound, the AF phase is 
observed at a much larger value of doping, when x ~ 0.03, where the spin-gap is already 
unobservable, (|||; whereas in the case of Si-doped compound, the AF phase is observed 
at y ~ 0.007, with co-existence of AF phase and spin-Peierls phase observed J/J. The much 
larger value of doping required in generating a AF state in the Zn-doped compound is con- 
sistent with the absence of inter-chain coupling between impurity- induced local moments in 
Zn-doped compound. The disappearance of spin-Peierls order in the AF phase also suggests 
that the AF phase may be formed by ordering of all spins in the Zn-doped compound, but 
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not ordering of impurity- induced local moments. Notice that co-existence of long-range an- 
tiferromagnetic and spin-Peierls order in the Ge-doped compound is a natural prediction of 
our theory if the AF state is formed by ordering of impurity-induced local moments, since 
the local magnetic moments can exist only in the presence of spin-Peierls order in our theory. 

It is also interesting to compare the low-energy properties of Si-doped spin-Peierls and 
the Zn-doped two-leg ladder compounds. In particular, interesting differences between the 
two compounds is expected at the intermediate temperature regime, when the intra-chain 
coupling is not important. In this case, the low energy physics of both compounds can be 
described by model of S — 1/2 Heisenberg spin chain with random nearest neighbor interac- 
tion, with one crucial difference: in the case of CuGei^ySiyO^, the effective nearest neighbor 
interaction has random magnitude but is always positive, whereas in the Sr{Cui- x Zn x )20^ 
compound, the effective nearest neighbor interaction has both random magnitude and sign 
pT3| . In the first case, the ground state of the system is described by a random singlet phase 



20| , |2~T|| which is characterized by formation of singlets between spins that may be far apart 



from each other. In the second case, effective spins with magnitude > 1/2 are generated 



at low energy because of ferromagnetic coupling [E3] and drives the system towards a fix 



point very different from the above [ 22"] . For example, it was found that the uniform mag- 
netic susceptibility x ~ T~ l in the random sign and magnitude model with different Curie 



constant at low and intermediate temperatures |T3| whereas \ ~ T l \lnT\ 2 in the random 



singlet phase. The specific heat C is also found to be very different in the two cases, with 
C ~ |ZnT|~ 3 in the random singlet phase and C ~ T 2a \lnT\ in the second case with a ~ 0.22 
2~2fl . Experimental tests of these differences are suggested. 

The interesting low energy behaviour associated with impurity-doped quasi- one dimen- 
sional spin-gap systems lead us to ask whether similar interesting phenomena can be ob- 
served in other quasi- ID spin systems. Indeed, it has been suggested that similar interesting 
physics can be observed in half-integer spin chains with S > 1 |17|| , which have bulk low 
energy physics identical to S = 1/2 spin chains It was predicted that local moments 



coupled to each other with effective interaction J' ~ l/L\lnL\ will form at ends of open 
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chains of length L. p7|j23|] . Following similar analysis as in Refs |13],[T7]], it can be shown 
that analogous low energy behaviour can be observed also in the three-leg Heisenberg ladder 
compound, upon substitution of Cu by Zn. It is convenient to view the three-leg ladder 
compound as three coupled antiferromagnetic spin chains. We find that replacement of Cu 
at the central chain of the three-leg ladder compound by Zn results at different physics 
from when the Cu at one of the side-chains is replaced p4| . First we consider replacing one 
single Cu atom by Zn. In the case when the Cu is located at one of the side-chains, we 
find that the Zn impurity effectively 'breaks' the ladder and the low energy physics of the 
system is the same as two uncoupled 5=1/2 spin chains; whereas in the case when the 
Cu atom is located at the central chain, we find that the Zn atom introduces effectively a 
5 = 1 impurity into the spin-ladder and the low energy physics is that of two S = 1/2 spin 
chains coupled antiferromagneticaly to the 5 = 1 impurity spin from both side. For finite 
concentration of Zn impurities, the low energy physics of the system can be described by 
5 = 1/2 spin chains coupled to randomly located 5 = 1 impurities (assuming that at least 
part of the Cu atoms at central chain are replaced by Zn) and the low energy behaviour of 
the system is determined by competition between Kondo screening of the impurity spins by 
the efffective 5=1/2 spin chain and antiferromagnetic ordering of impurity- induced 
5 = 1 local moments, which may lead to very different low temperature behaviour compared 
with the spin-Peierls or two-leg ladder systems ||24)| . 

I acknowledge helpful discussions with H. Fukuyama, A. Azuma, N. Nagaosa, A. Furusaki 
and M. Sigrist. I also thank the Yukawa Institute, Kyoto, where part of this work was done. 
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